home *** CD-ROM | disk | FTP | other *** search
/ io Programmo 60 / IOPROG_60.ISO / soft / c++ / gsl-1.1.1-setup.exe / {app} / src / cblas / test_her.c < prev    next >
Encoding:
C/C++ Source or Header  |  2002-04-18  |  4.4 KB  |  180 lines
  1. #include <gsl/gsl_test.h>
  2. #include <gsl/gsl_ieee_utils.h>
  3. #include <gsl/gsl_math.h>
  4. #include <gsl/gsl_cblas.h>
  5.  
  6. #include "tests.h"
  7.  
  8. void
  9. test_her (void) {
  10. const double flteps = 1e-4, dbleps = 1e-6;
  11.   {
  12.    int order = 101;
  13.    int uplo = 121;
  14.    int N = 1;
  15.    int lda = 1;
  16.    float alpha = 1.0f;
  17.    float A[] = { 0.188f, 0.856f };
  18.    float X[] = { -0.832f, -0.151f };
  19.    int incX = -1;
  20.    float A_expected[] = { 0.903025f, 0.0f };
  21.    cblas_cher(order, uplo, N, alpha, X, incX, A, lda);
  22.    {
  23.      int i;
  24.      for (i = 0; i < 1; i++) {
  25.        gsl_test_rel(A[2*i], A_expected[2*i], flteps, "cher(case 1410) real");
  26.        gsl_test_rel(A[2*i+1], A_expected[2*i+1], flteps, "cher(case 1410) imag");
  27.      };
  28.    };
  29.   };
  30.  
  31.  
  32.   {
  33.    int order = 101;
  34.    int uplo = 122;
  35.    int N = 1;
  36.    int lda = 1;
  37.    float alpha = 1.0f;
  38.    float A[] = { 0.188f, 0.856f };
  39.    float X[] = { -0.832f, -0.151f };
  40.    int incX = -1;
  41.    float A_expected[] = { 0.903025f, 0.0f };
  42.    cblas_cher(order, uplo, N, alpha, X, incX, A, lda);
  43.    {
  44.      int i;
  45.      for (i = 0; i < 1; i++) {
  46.        gsl_test_rel(A[2*i], A_expected[2*i], flteps, "cher(case 1411) real");
  47.        gsl_test_rel(A[2*i+1], A_expected[2*i+1], flteps, "cher(case 1411) imag");
  48.      };
  49.    };
  50.   };
  51.  
  52.  
  53.   {
  54.    int order = 102;
  55.    int uplo = 121;
  56.    int N = 1;
  57.    int lda = 1;
  58.    float alpha = 1.0f;
  59.    float A[] = { 0.188f, 0.856f };
  60.    float X[] = { -0.832f, -0.151f };
  61.    int incX = -1;
  62.    float A_expected[] = { 0.903025f, 0.0f };
  63.    cblas_cher(order, uplo, N, alpha, X, incX, A, lda);
  64.    {
  65.      int i;
  66.      for (i = 0; i < 1; i++) {
  67.        gsl_test_rel(A[2*i], A_expected[2*i], flteps, "cher(case 1412) real");
  68.        gsl_test_rel(A[2*i+1], A_expected[2*i+1], flteps, "cher(case 1412) imag");
  69.      };
  70.    };
  71.   };
  72.  
  73.  
  74.   {
  75.    int order = 102;
  76.    int uplo = 122;
  77.    int N = 1;
  78.    int lda = 1;
  79.    float alpha = 1.0f;
  80.    float A[] = { 0.188f, 0.856f };
  81.    float X[] = { -0.832f, -0.151f };
  82.    int incX = -1;
  83.    float A_expected[] = { 0.903025f, 0.0f };
  84.    cblas_cher(order, uplo, N, alpha, X, incX, A, lda);
  85.    {
  86.      int i;
  87.      for (i = 0; i < 1; i++) {
  88.        gsl_test_rel(A[2*i], A_expected[2*i], flteps, "cher(case 1413) real");
  89.        gsl_test_rel(A[2*i+1], A_expected[2*i+1], flteps, "cher(case 1413) imag");
  90.      };
  91.    };
  92.   };
  93.  
  94.  
  95.   {
  96.    int order = 101;
  97.    int uplo = 121;
  98.    int N = 1;
  99.    int lda = 1;
  100.    double alpha = 0.1;
  101.    double A[] = { 0.257, 0.326 };
  102.    double X[] = { 0.319, -0.009 };
  103.    int incX = -1;
  104.    double A_expected[] = { 0.2671842, 0.0 };
  105.    cblas_zher(order, uplo, N, alpha, X, incX, A, lda);
  106.    {
  107.      int i;
  108.      for (i = 0; i < 1; i++) {
  109.        gsl_test_rel(A[2*i], A_expected[2*i], dbleps, "zher(case 1414) real");
  110.        gsl_test_rel(A[2*i+1], A_expected[2*i+1], dbleps, "zher(case 1414) imag");
  111.      };
  112.    };
  113.   };
  114.  
  115.  
  116.   {
  117.    int order = 101;
  118.    int uplo = 122;
  119.    int N = 1;
  120.    int lda = 1;
  121.    double alpha = 0.1;
  122.    double A[] = { 0.257, 0.326 };
  123.    double X[] = { 0.319, -0.009 };
  124.    int incX = -1;
  125.    double A_expected[] = { 0.2671842, 0.0 };
  126.    cblas_zher(order, uplo, N, alpha, X, incX, A, lda);
  127.    {
  128.      int i;
  129.      for (i = 0; i < 1; i++) {
  130.        gsl_test_rel(A[2*i], A_expected[2*i], dbleps, "zher(case 1415) real");
  131.        gsl_test_rel(A[2*i+1], A_expected[2*i+1], dbleps, "zher(case 1415) imag");
  132.      };
  133.    };
  134.   };
  135.  
  136.  
  137.   {
  138.    int order = 102;
  139.    int uplo = 121;
  140.    int N = 1;
  141.    int lda = 1;
  142.    double alpha = 0.1;
  143.    double A[] = { 0.257, 0.326 };
  144.    double X[] = { 0.319, -0.009 };
  145.    int incX = -1;
  146.    double A_expected[] = { 0.2671842, 0.0 };
  147.    cblas_zher(order, uplo, N, alpha, X, incX, A, lda);
  148.    {
  149.      int i;
  150.      for (i = 0; i < 1; i++) {
  151.        gsl_test_rel(A[2*i], A_expected[2*i], dbleps, "zher(case 1416) real");
  152.        gsl_test_rel(A[2*i+1], A_expected[2*i+1], dbleps, "zher(case 1416) imag");
  153.      };
  154.    };
  155.   };
  156.  
  157.  
  158.   {
  159.    int order = 102;
  160.    int uplo = 122;
  161.    int N = 1;
  162.    int lda = 1;
  163.    double alpha = 0.1;
  164.    double A[] = { 0.257, 0.326 };
  165.    double X[] = { 0.319, -0.009 };
  166.    int incX = -1;
  167.    double A_expected[] = { 0.2671842, 0.0 };
  168.    cblas_zher(order, uplo, N, alpha, X, incX, A, lda);
  169.    {
  170.      int i;
  171.      for (i = 0; i < 1; i++) {
  172.        gsl_test_rel(A[2*i], A_expected[2*i], dbleps, "zher(case 1417) real");
  173.        gsl_test_rel(A[2*i+1], A_expected[2*i+1], dbleps, "zher(case 1417) imag");
  174.      };
  175.    };
  176.   };
  177.  
  178.  
  179. }
  180.